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We present a calculation of the B and D meson decay constants in lattice QCD with two 
(Nf — 2) flavours of light dynamical quarks, using an 0(a)-improved Wilson action for both light 
and heavy quarks and a renormalization-group improved gauge action. Simulations are made at 
three values of lattice spacing a = 0.22, O.f 6, O.lf fm and four values of sea quark mass in the range 
mps/mv ~ 0.8-0.6. Our estimate for the continuum values of the decay constants are fsd = 
208(10) (ll)MeV, / Ss = 250(10) (13) (± e )MeV, f Dd = 225(14) (14)MeV, f D . = 267(13) (17) (±J°)MeV 
for Nf — 2 where the statistical and systematic errors are separately listed, and the third error for 
/b„ and fos show uncertainty of determination of strange quark mass. We also carry out a set of 
quenched simulations using the same action to make a direct examination of sea quark effects. Taking 
the ratio of results for Nf — 2 and Nf = 0, we obtain f B ^ 2 /f B ^° = 1.11(6), fgf 2 //b 3 /_ ° = 
1.14(5), fo f d =2 /fn f 4 = ° = 1-03(6), /ff =2 //^ff = ° = 1.07(5). They show a 10-15% increase in the 
Nf — 2 results over those of Nf = for the B meson decay constants, while evidence for such a 
trend is statistically less clear for the D meson decay constants. 

PACS: 12.15.Hh, 12.38. Gc, 13.20.-v 



I. INTRODUCTION 



The accurate determination of the CKM matrix elements is one of the most important tasks of flavour physics. The 
Standard Model prediction of the unitarity of the matrix still has to be tested, especially for the unitarity relation 
involving the most off-diagonal elements, which contain the source of the CP violation in the Standard Model. 

Two of the matrix elements, \Vtd\ and |Vt s |, can be extracted from experimental data of the oscillation frequency 
Am, of B q — B q systems (q denotes either d or s quark) through the relation Q| 

Am, = ^M^So(x t )rj B M B J 2 Bq B Bq \V tq V t l\ 2 , (1) 

where the factors other than fg q B Bq are known either experimentally or through perturbative calculations in QCD. 
The non-perturbative coefficients f Bq and B Bq are defined as 

(V\hul?>q\B q (p))=if Bq p a , (2) 

and 

6 p / \ ( B g\hn(i - 75)g^( 1 ~ i5)q\B q ) 

B Bq = RbW 8 f 2 M 2 ' W 

3 ^ B q B q 



*Present address: Department of Physics, Columbia University, 538 West 120 th St., New York, NY 10027, U.S.A. 
^Present address: Department of Biochemistry and Molecular Biology, University College London, Darwin Building, Cower 
St., London, WC1E 6BT, England, U.K. 



1 



where denotes a renormalization group factor to eliminate the variation due to the scale /x where the four-quark 

operator 67^(1 — 75)967^(1 — 75)9 is defined. In this paper we shall focus on the decay constants /b , leaving the bag 
parameter Bs q for future studies. 

Experimentally, the Cabibbo-allowed leptonic decay D s — > tv t has been measured and the recent values for fjj g are 
285±20±40 MeV (ALEPH @) and 280±19±44 MeV (CLEO §). On the other hand, a measurement of the decay 
constant fs is difficult, since B + — ► l + Di is Cabibbo-suppressed in the Standard Model. Hence fs q has to be provided 
from theory, while fu can be used to check the calculational method. 

The calculation of these decay constants have been carried out extensively in the quenched approximation in lattice 
QCD, where vacuum polarization effects are neglected in order to reduce the computational requirements. A recent 
summary of these attempts is given in Refs. Although the approximation provides a useful first step in a lattice 

QCD determination of the decay constants, the size of the resulting systematic error is not clear. In the quark potential 
model, the decay constant is proportional to the wave function at the origin squared and the potential at such short 
distance scales can be expressed in terms of the running coupling constant. Therefore, one can demonstrate, on a 
heuristic level, that the decay constant is affected by the number of flavours Nf, where Nf — corresponds to the 
quenched approximation. Additionally, a recent study of the light hadron spectrum in the quenched approximation 
indicates a deviation of about 10% from experiment |6j. For the decay constants, quenched chiral perturbation theory 
suggests that the deviation introduced by the approximation may be significant. 

The elimination of this approximation is numerically highly intensive and has become realistic only recently. The 
MILC collaboration Q and Collins el al. jl0| have performed first calculations of the decay constants on the lattice 
with two degenerate sea quark flavours, and found an indication that fs is considerably larger in the presence of sea 
quarks. In these studies, however, the discretization of the sea quarks is defined using the staggered fermion action, 
which is different from that used for the light valence quark (Wilson fermion in Ref. |)| and the O(o)-improved (clover) 
fermion in Ref. It would complicate the discretization error in the results. In fact, the observed a, (the lattice 

spacing), dependence in Jb found in Ref. ^| is rather different between quenched and unquenched calculations, even 
though the formulations for valence heavy and light quarks are the same. 

In our work we apply a consistent formulation where for both sea and light valence quarks we use the same action, 
and study the a dependence by performing three sets of two-flavour calculations at a ~ 0.22, 0.16 and 0.11 fm. For 
comparison, we carry out quenched calculations at ten different values of a covering the range studied in the two 
flavour calculations. We employ the 0(a)-improved quark action |ll]] for both sea and valence light quarks. The 
same action is used for the heavy quark, applying the non-relativistic reinterpretation of Ref. p"2| . The gauge field 
is described by a renormalization group improved action [T^| , which reduces the descretization error on the coarse 
lattices on which our calculations are made. 

The rest of this paper is outlined as follows. In Sec. || we discuss the lattice actions and the formulation to 
treat heavy quarks. The computational details involved in the calculation are described in Sec. Ill, and our analysis 
procedures in Sec. IV. We present the results in Sec. [y] where we discuss in particular how we estimate the values 
in the continuum limit and their errors, and make a comparison between the Nf — and Nf — 2 results. Our 



conclusions are summarized in Sec. VI 



II. LATTICE ACTIONS 



A. Light sector 

The renormalization group (RG) improved gauge action we employ takes the form [ fl3f 

S f ^(coE^xi+^E 1 ^)' (4) 

where W± x i and W± X 2 are the Wilson loops of size lxl and 1x2 respectively, and the sums run over all possible 
sites and orientations. The parameter (3 is related to the bare gauge coupling g% through (3 = 6/(70 ■ The coefficients 
Co and c\ are defined as 

c = 3.648, (5) 
ci = i (1 - c ) = -0.331, (6) 

which are chosen so as to approximate the renormalization gro up trajectory in two dimensional operator space. 
For quarks we employ the 0(a)-improved (clover) action |ll[| defined by 
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1>v, (?) 
where is the standard Wilson formulation of the Dirac fermion matrix 

= S xy - K^2{(1 - 7/1 )C4,A+A,a + (1 + iMjw+fi} (8) 
and the matrix F^ v is the simplest definition of the field strength, 

Fpu = ^7 (//«/ — /^l/) ? (9) 

where / M „ is the standard clover-shaped definition of the gauge field strength. The leading discretization error in the 
Wilson fermion action (csw = 0) is removed by appropriately tuning the parameter csw- We apply a mean field 
approximation c$w = P~ 3 ^ 4 , where P—(Wi x i), and perturbative expansion at one-loop P = 1 — 0.1402<?q is used 
to evaluate P. With this choice, the leading contributions among remaining discretization errors are 0(a s a) and 
0(a 2 ) for light quarks. There is reasonable agreement between P measured on the lattice and the above perturbative 
definition with the difference being at worst 8% p^j. Furthermore, there is also good agreement between the above 
definition of csw and the one-loop computed value |Q. 

The efficacy of this choice of actions over the standard action has been demonstrated in Ref. |H by examining 
the rotational invariance of the static potential and the scaling behaviour of the light hadron spectrum. In using the 
clover fermion action, we also note that care must be taken in defining currents, which will be discussed below. 
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B. Heavy Quarks 



It seems implausible to examine hadrons containing heavy quarks with mass toqo > 1 on a lattice with the spacing 
a, as one expects the discretization effects to become uncontrollably large for such large masses. However, this is not 
necessarily true for heavy-light mesons. The spatial momentum of the light degrees of freedom in the heavy-light 
system is controlled by the QCD scale Aqcd rather than the much larger heavy quark mass scale. In the limit of 
infinite mo, the heavy quark mass decouples from the dynamics of the system, and the Heavy Quark Effective Theory 
(HQET) fllq ] becomes a good approximation. At a finite rriQ, the correction may be incorporated as an expansion in 
1/rriQ, which is a basis of the non-relativistic QCD (NRQCD). 

On the lattice, it is straightforward to formulate the static jlTj and NRQCD jl8|,[l9) actions, and a number of 
(quenched) calculations of fg have been performed using them. Another formulation to realize the idea of HQET on 
the lattice |t2| is also useful, as it uses the same relativistic form of the quark action as that for light quarks except 
that the bare heavy quark mass mo may be taken to be arbitrarily large. 

For the heavy-light system, where the typical spatial momentum of the heavy quark is small compared to the inverse 
lattice spacing, one can construct an effective Hamiltonian starting from a relativistic action, 



H w * 



Mi + 70 A) 



D 2 



/SB 



2M 2 2Mb 



7o- 



[7 -D.7-E] 
8M E 



(10) 



where D is the covariant derivative, X the Pauli spin matrices, and B and E are the chromomagentic and chro- 
moelectric fields respectively and an expansion in small spatial momentum or equivalently in aD on the lattice is 
performed. This Hamiltonian is equivalent to the standard non-relativistic Hamiltonian if the "mass" parameters Mi, 
M%, Mb and Me are equal to each other. Those are, however, different functions of am and not necessarily equal 
to each other, unless the parameters in the initial relativistic action are appropriately tuned. The strategy suggested 
in Ref. |0| is, therefore, to take the action as an effective theory to generate the dynamics described by ([lO|). The 
appropriate mass parameter in the non-relativistic effective theory is the "kinetic" mass M2, while the "pole" mass 
Mi does not affect the dynamics of heavy quark and plays merely a role of energy shift in this formalism. To obtain 
a correct action at order 1/M, the mass parameter which characterizes the spin-magnetic interaction Mb must be 
equal to M2, which is satisfied for the 0(a)-improved (clover) action up to perturbative corrections. On the other 
hand, there is no tunable parameter in the clover action to make Me equal to Mi and Mb, so that the contributions 
of 0(1/M 2 ) and higher are not correctly described by the clover action. 
At tree level, the kinetic mass Mi of the heavy quark is given by |L2] 
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aM 2 = ( 77^ r + 1 , (11) 

\amo(2 + arrio) 1 + amo J 

where the bare mass amo is defined as amo — — The one- loop relation is also known [2p|-|2^1 for the standard 
plaquette gauge action but not for our choice of the action. Hence we employ the tadpole improvement [p3| of the 
above relation, which is obtained by simply replacing amo with 8K c amo. 

The heavy-light meson mass aMj^QET defined in the HQET is then obtained as ]24| , |2l| 

vMMqet = aM™ le + (aM 2 Q - aM«), (12) 

from the pole mass aM po i e extracted from the exponential fall off of the heavy-light propagator. The superscript M 
or Q in (Q) distinguishes the mass of the heavy- light meson (M) from the heavy quark mass (Q). The parameters 
aM® and aM® are the tree-level defined pole and kinetic masses of the heavy quark. 

An alternative way to obtain the heavy-light meson mass is to measure its energy-momentum dispersion relation 
and fit with the form E(p) — M po i e + p 2 / (2Mki n ) + 0(p 4 ) to extract the "kinetic" mass Mkin (as employed in pcf). 
Unfortunately, for the lattices which were used to quote our final results the statistical ensemble was not large enough 
to obtain an accurate measurement of M/-j n . For this reason, this choice of the kinetic mass will not be further 
discussed here. 

The axial current to be measured should also be modified to obtain the results correct at 0(1/M) according to 

h -> (1 - adi7 • D)/i, (13) 

where h is the heavy quark field and equivalently for h, and the parameter d\ is a function of amo- At the tree level, 
it is given by p2[ 

a«i = j- r — -—rr, (14) 

am (2 + am ) 2aM 2 

and the axial vector current for heavy-light mesons, correct to 0(1/M), takes the form 

A^x) = /(x)757 P /i(x) - adiF(x)757 M 7 • Ah(x). (15) 

where I is the light quark field. The tadpole improvement of d\ may be applied again with the replacement amo ~~ * 
8K c amo- 

The following point should also be noted. The action being used is still a relativistic action and as the lattice spacing 
becomes smaller, it is expected that theory should smoothly cross over to a fully relativistic theory. That means the 
mass parameters M's become identical as amo decreases. The lattice spacing dependence of physical quantities, such 
as fs and /d, is, however, highly nontrivial unless mo is much smaller than 1/a, and the continuum extrapolation 
in such a situation would not be justified with any simple ansatz, e.g. linear or quadratic in a. The formulation is, 
therefore, treated as an effective theory (like NRQCD), and the discretization error should be reasonably small at 
fixed a in order to obtain reliable results. 

Despite the caveat of the preceding paragraph, this approach has been successfully implemented in the quenched 
approximation in Refs. [p5| , p6| using the plaquette gauge action. Since we use a gauge action which has been unused 
in the previous heavy quark calculations it is important for us that we repeat the calculation in the quenched 
approximation in order to see if the quenched results obtained with the "standard" plaquette gauge action are 
reproduced. 



III. COMPUTATIONAL DETAILS 



A. Gauge Fields 

Gauge configurations were generated for Nf—0 and Nf=2 using the RG improved gauge action and the 0(a)- 
improved Wilson quark action as discussed in Sec. II A. Technical details on the configuration generation for Nf = 2, 
carried out with the Hybrid Monte Carlo algorithm, are described in our dynamical QCD calculations papers [p7|-p9[ . 

In the Nf—2 calculations, we performed three sets of calculations at bare gauge couplings j3 = 1.8, 1.95 and 2.1, 
which correspond to the lattice spacing a ~ 0.22, 0.16 and 0.11 fm respectively. The lattice size is 12 3 x24 (/3=1.8), 
16 3 x32 (1.95) and 24 3 x48 (2.1), with which the physical volume is approximately (2.5 fm) 3 . For each set, we carried 
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out runs at four values of sea quark mass in order to take the chiral limit of sea quark. The four sea quark masses 
are tuned so that the pseudoscalar-to-vector mass ratio mps/mv becomes roughly 0.80, 0.75, 0.70 and 0.60, which 
correspond to the range of quark mass of 3 ~ 0.5 times physical strange quark mass. The simulation parameters are 
listed in Tabic |, where the number of HMC trajectories is also shown. We note that at (3 = 2.1 the configurations 
analyzed constitute the first half of the ensemble for each sea quark mass. The full set of configurations is used at 
j3 = 1.95 and 1.8. The measurements are performed on configurations separated by 10 HMC trajectories at /3 = 1.8 
and 1.95 and by 5 trajectories at (3—2.1. The statistical analysis is done using the jackknife method in order to take 
the correlation of successive trajectories into account. The bin size is 50 trajectories for all Nt = 2 runs, which has 
been determined to be a suitable length for eliminating autocorrelations Q . 

The lattice spacing is determined for each (3 value using the p meson mass as input at the physical sea quark limit. 
The chiral extrapolation of light hadrons is discussed in |2l] 2j| . The lattice spacings are listed in Table ||. 

In order to see the sea quark effect consistently using our choice of gauge and quark actions, we prepared ten sets of 
the quenched (Nf=0) gauge configurations. The values of (3 are chosen so that the string tension matches with each 
of full QCD configurations at (3 = 1.95 or 2.1 at four sea quark masses and also in the chiral limit. For calculating 
lattice spacing and hence the physical value of the decay constants, the p meson mass is used as input in conjunction 
with the vector masses measured on the lattice extrapolated to the light quark masses. The detail of our parameter 



choice in the quenched runs is summarized in Table HI 



B. Valence quarks 

The heavy and light quark propagators arc calculated on each set of the gauge configurations for the O(o)-improved 
Wilson action with the same choice of csw as used in the configuration generation. For each set of gauge configurations, 
eight values of the heavy quark mass are chosen so that their HQET mass ([l2|) lie roughly on the interval of the b 
and c quark masses. 

The light quark mass on the dynamical configurations is the same as their sea quark mass. In addition, we choose 
another quark mass for each set of configurations so that it satisfies mpsjvny — 0.688. To compute any of the 
observables at the strange quark mass, the relevant observables are interpolated to the strange quark mass defined 
from the mass of the K or tp. The light quark masses for Nf — are chosen to take values approximately the same 
as those for the equivalent lattices (i.e. those lattices with matched string tension) for Nf = 2. 

The gauge configurations are fixed to the Coulomb gauge with a global maximum residue for Tr(diAi) 2 set to 10~ 14 
or less. The light quark propagators are solved with local sources while the heavy quark propagators are computed 
with local and smeared sources. The smearing is made with the exponential function ^4exp(— Br), with the mean 
radius l/B chosen to approximately reproduce the heavy-light wave function. The parameters A and B are listed in 



Table [TV 

Both light and heavy quark propagators are obtained with a solver based on the BiCGStab algorithm. For large 
values of heavy quark mass, stopping the solver if the residue becomes smaller than some minimum is not sufficient 
for obtaining the solution at large time separations. In this case, the iteration of the solver is applied a minimum of 
2 x T times, where T is the temporal extent of the lattice, before applying the maximum residue criterion. 

C. Heavy-light current 

We compute the correlation functions constructed from the following operators 

P(x) = l{xh 5 h{x), (16) 
A(x) = l(xh 5l0 h(x), (17) 
8A{x) =Z» 7 s7o7- AA(a). (18) 



The heavy and light quark fields h and I are normalized with y 1 — 3K/4K C , which is motivated with the nonrelativistic 
interpretation |L2j together with the tadpole improvement |2^]. The derivative current 5 A is used to construct the 
modified current according to ([H]), and A is the discretised covariant derivative defined as 

Aih(x) = ^ Ui(x)h{x + i)-U}(x-i)h{x-i) . (19) 
Specifically, we measure the correlation functions 
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EAM^mm, Eg(M^) p m), (20) 
£ 2 <M(af,t)pt(o)>, E s (tM*,t)Plm, 

where the subscripts S and L on the pseudoscalar operators indicate whether smeared or local operators are employed. 
The axial current A and the derivative current are always local. 

IV. ANALYSIS 
A. Correlators 

The correlation functions defined above take the following form for large Euclidean time separation (we take a = 1 
for simplicity) 

J2(PL&t)Pl LS) (0)) = ZFL 2m L ' S) e- M *° l ° T/2 cosh (M pole (T/2 - t)) 

X 

+ Zk 2w' S) *- MUT/ * cosh «ol*(T/2 ()), (21) 
J2(A(x,t)P^ s) (0)) = ZA \ P ^ e~ M ^ 2 sinh (M pole (T/2 t)) 

X 

+ e ~ K ° leT/2 Sinh ( M Pole(T/2 - t)), (22) 

^2(SA(x,t)Pl LjS) (0)) = A ^ S) e - M >°" T ' 2 sinh (M pole (T/2 - t)) 

>M ^El e -M' poU T/2 sinh {M ' pole (T/2 - t)), (23) 

where T is the temporal extent of the lattice, and M and M' are masses of the ground and the first excited pseudoscalar 
states respectively. The masses extracted from the t dependence of the correlation functions are the pole masses, while 
the M's appearing in the denominator come from the normalization of states, and their definition need not be specified 
for calculating the combination of fpyM. 
The matrix elements Z are defined as 

Z PlLS) = (0|P (i)S) (0)|P(0)), (24) 
Z A = (0|A(0)|P(0)>, (25) 
SZ A = (0|cL4(0)|P(O)), (26) 

where |P(0)) represents the heavy-light pseudoscalar meson state at rest. The primed quantities are defined in a 
similar manner for the first excited state. 

We carry out a simultaneous fit of the three correlators ([2l|)-(]23]). Formally, for large Euclidean times, the con- 
tribution of the excited state will be negligible. However, it is included in the fit so as to use a wider range of 
Euclidean times and reduce the size of the statistical error. For those sets of configurations with lattice volumes of 
size 12 3 x 24 and 16 3 x32 the sample size is large enough to perform a correlated fit of the local and smeared source 
data, where the correlation among different time slices are also taken into account. For the largest lattice 24 3 x48, 
such a correlation matrix appeared to be too large to achieve a stable fit with our statistics. We, therefore, use the 
uncorrelated fit throughout our statistical analysis, and check that the results are unchanged within statistical errors 
with the correlated fit when it is possible. 

The fit criteria we apply for selecting the fit range are as follows |3C]]: (i) The quality of fit, Q, should be acceptable, 
e.g., Q > 0.1. (ii) The results for the chosen fit range should agree to within one standard deviation of the results 
when the minimum time slice is increased or decreased by one time-slice, (iii) There should be agreement between the 
ground state results obtained using a single-exponential fit and a double-exponential fit. This condition increases our 
confidence that higher state contamination is eliminated, (iv) In the double-exponential fit the ground and excited 
state energies must be statistically resolvable, i.e., there must be more than one standard deviation between their 
central values (since we expect the physical states to be distinctly seperated). 
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The effective mass plots for the (PP) correlators, together with fit curves, are shown for a typical heavy-light meson 
mass in Figures [j]-|3] (Nf = 2 case) and in Figure |]( quenched case). 

B. Heavy-light decay constant 

The heavy-light decay constant fp is obtained through 

a 3 / 2 (fpVM) = Z A -^=(Z A -ad 1 8Z A ). (27) 

In the massless limit, the renormalization constant Z A was previously calculated perturbatively to one-loop order for 
the RG- improved action fjlf . Here we use a recent extension of this result to finite heavy quark masses made by 
K.-I. Ishikawa et al. |32]. The results can be expressed in the form 



Z A = 1 



p A log(am ) 

TT 



(28) 



where the one-loop coefficient po = p A — ^ log(amo) is plotted in Figure || as a function of amo, for the cases when 
d± takes the tree-level value and when it is ignored. 

It is well known that perturbative expansions in lattice QCD are ill-behaved when one uses the bare coupling con- 
stant i?g , and the use of some renormalized coupling defined through short-distance quantities gives a more convergent 
expansion p3| . Since a two- loop calculation of short-distance quantities necessary to define an appropriate renormal- 
ized coupling is not yet available for the RG gauge action, we use the continuum MS coupling as an alternative. 

The one-loop perturbative relation between the bare and MS couplings for the RG improved gauge action and the 
0(a)-improved Wilson quark action is known as J3l| 

1 § + 0.1000 + 0.0315^/. (29) 



gf— (u = l/a) 6 



The tadpole improvement p3fl may be applied to reduce the ultraviolet dominated pieces from the perturbative 
expansions by reorganizing the above relation as 

' (c P -8c x R)~- 0.1006 + 0.0315iV/, (30) 



where P = (VFixi) and R = (Wi X 2) are the expectation values of plaquette and 1x2 rectangle jlj|, and the one-loop 
expressions P = 1 — 0.1402<? 2 and R = 1 — 0.2689g 2 are used to obtain the modified one-loop coefficient in (|30|). The 
values of gj—(p = 1/a) obtained with this formula are 3.162, 2.812 and 2.562 at (3—1.8, 1.95 and 2.1 respectively for 



the Nf = 2 lattices. The same quantity for the quenched lattices is listed in Table [II. 

In Figure ^| we plot Z A as a function of the bare heavy quark mass for the plaquette and RG-improved actions for 
an inverse lattice spacing of around 1.8 GeV (/3=5.9 for the Wilson and P— 2.528 for the RG action in the quenched 
approximation). In contrast to the large one-loop correction of order —20% for the case of the plaquette gauge action, 
Z A is close to unity for the RG-improved action. 

V. RESULTS 
A. Effect of field rotation to the heavy-light current 



We first examine the effect of field rotation (13), which is reflected in ( |15[ ) and ( |27| ) as a correction proportional to 
d\. An order counting suggests the size of the correction of the order of adiAQCD, which is about d±xl5% at 1/a w 2 
GeV if Aqcd — 300 MeV. Since the tree-level coefficient d\ given by (|lj) is smaller than 0.1 for any value of the bare 
quark mass amo, the size of the correction is naively estimated to be 0(2%). 

In Figure @ we plot the quantity /rotated^/ junrotated _ i as a function of the meson mass for Nf = 2 and Nf = 0, 
where ^rotated j ncmc j e s the rotation term while it is ignored in jumotated Q are mus t be exercised in this comparison 
to use the appropriate renormalization factors Z A for the rotated and unrotated currents shown in Figure ^ since the 
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diagram originating from the rotation term should be excluded for calculating junrotatea_ rp^e i^ice spacing for 
(3 = 2.575 at Nf = is approximately equal to the lattice spacing, extrapolated to the chiral limit, for j3 — 2.1 at 
Nf — 2, which allows a more relevant comparison of the ratios. As one can see, the magnitude of correction is of the 
order of 3-7%, which is larger than our expectation and cannot be ignored. 

The large magnitude of this correction may partly originate from a power divergence of the matrix element of the 
higher dimensional operator 5 A defined by (Cq), with which the naive order counting of 0(o,Aqcd) is changed to a size 
of O(l). In principle this power divergence should be compensated by that in the perturbative matching. However, 
at the one-loop order in the calculation of Ref. [B2|, the compensation is incomplete. 



B. Extrapolation to physical quark masses 

In order to obtain the heavy-light decay constant fp^/M for the physical mass of B( s ) and Z?( s ) mesons, we fit the 
data with the following form 

a 3 / 2 $ P = A + A ia m q + A 2 (am q ) 2 + — [B + B ia m q ] + T^y C , (31) 

where we define the renormalization group invariant decay constant 4>p as 

a 3/ 2$p = f^iMl) 27/30 a^(f P VM), (32) 
V a s [Mb ) / 

with /?o = 11 — f-W/. The light quark mass is defined as am q = i — where K c denotes the value at which 

pion mass made of sea quarks vanishes, and the HQET mass definition ([l2]) is used for the heavy-light meson mass 
M . The renormalization group factor (a s (M) / a s (MB)) 2 ^° is evaluated with a two-loop running coupling coefficient 
adopting Aqcd = 300 MeV for both Nf = 2 and Nf = 0. We have checked that the resulting decay constants are 
stable well within statistical errors under a variation of Aqcd by a factor two. 



The form (31) is a truncated expansion of the matrix element in 1/aM and in am q . It is possible to include higher 
order terms; however, the resulting fit coefficients are statistically not well determined, and we do not include such 
terms in our analyses. 

In determining fg d in the Nf — 2 case, we only employ the matrix elements where the sea and valence light 
quark masses are matched. For fg s we interpolate, at each sea quark mass, the matrix element in the valence light 
quark mass to the physical strange quark determined using the partially quenched analysis p7|-p9|. The values of 
the hopping parameter K s corresponding to the strange quark are listed in Tables || and III , for the K and (f> meson 



masses as physical input. The critical hopping parameter K c necessary for evaluating the light quark mass m q is also 
listed in these tables. 

The quenched data are analyzed with the same fit ansatz except that the term Ai is set to zero, as the number of 
light quark masses in this case precluded a quadratic fit. For fs B and fn„ the terms A\ and B\ are also set to zero, 
since there is no remaining light quark mass dependence once the strange quark mass is fixed. 

Fits with the form ( |3l|) are represented in Figures |[ |9| and [l0| for three lattice spacing for Nf = 2 and in Figure |ll| 
for Nf = 0. Data points and fit curves are plotted as a function of 1/aM for fixed am q , from which one can see that 
the ansatz ( |3l| ) represents the data quite well, except for a few points at /3=1.8. 

It is also illustrative to plot the data at fixed heavy quark masses as a function of am q , which is shown in Figures |l2] - 
[Ti] for Nf — 2 and Figure |l5| for Nf — 0. Since the results are given for fixed Kh, the heavy hopping parameter, 
we interpolated the curves for aM as a function of Kh and am q and hence reexpressed the coefficients of ([5l|) as a 
function of Kh and am q . For Nf = 2 we find clear curvature, which motivated us to introduce the term ^(am,) 2 in 
(|3l|). We, then, find good agreement of the fits to the data points. The fit parameters Ai, Bi and Co are summarized 
in Tables 0-|x| for each set of configurations. 

The B and D meson decay constants in physical units are obtained from (131]) with their physical meson masses as 



input, and the numbers are summarized in Tables XI and XII respectively. The lattice scale is set using the p meson 
mass. 



C. Discretization effect 

The decay constants are plotted as a function of a in Figures |iq(/B rf ), 0(/b s ), Q/dJ and |19|(/dJ. For Jb b 
(Figure |l7|) and fu a (Figure |l9|), we use the mass of K to define the strange quark mass (as a short hand, we will 
refer to this as fs s (K) and fo s {K)). 
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For the quenched data (Nf — 0), where ten data points are available, we observe a linear decrease as the lattice 
spacing decreases from awl GeV -1 to a w 0.8 GeV -1 , followed by an almost constant behavior within statistical 
fluctuations below a ~ 0.8 GeV -1 . We therefore fit the five data points for a< 0.8 GeV -1 as shown in the figures, and 
take this as our best estimate for the quenched value of the decay constant. 

There is no evidence that Nf = 2 data become independent of the lattice spacing. So we are not able to safely 
estimate the B meson decay constant from our data. We may discuss, however, our results in the following way. The 
slope of the decrease of Jb from (3 — 1.8 to 1.95 quite resembles that for Nf — for stronger couplings, while the 
decrease from (3 = 1.95 to 2.1 is somewhat reduced. If we suppose that the Nf — 2 data behave in a way similar to 
those for Nf — 0, the Nf — 2 data would have already reached the asymptotic flattening at around a » 0.7GeV _1 
and the data at (3 = 2.1 may be taken as a continue value. Since we cannot do better with the present data, we 
provisionally take the point at [3 = 2.1 as the continuum value, allowing for the possibility that the true value may 
be somewhat smaller than our estimate. From the shape of the beta dependence of the Nt — 2 data and their error 
bars, however, we may still safely conclude that the continum value does not undershoot that of Nf = 0. 

The extraction of the continuum limit for the D meson decay constant is more subtle, since we see a larger drop 
from j3 = 1.95 to (3 = 2.1 rather than (3 — 1.8 to 1.95. While we take the data at j3 = 2.1 as our provisional estimate 
for fjj, there is a possibility that the true value is smaller, although we expect that it is equal to or larger than 
quenched result from the experience with the B decay constant. Our safe conclusion is that the dynamical effect for 
D mesons is apprecaibly smaller than that for B mesons. 

We employ the same strategy as above for estimating the ratios fs s /fBd an d fDs/fod as shown in Figures ^ and 



D. Systematic errors 



We now examine the issue of systemtic errors in our results for the decay constants. For this purpose we list the 
possible leading order errors and estimate their magnitude using naive power counting. 

Generically these errors appear in three forms. First, we use tree-level mean-field estimates of the coefficients in 
the actions and currents and hence there will be radiative corrections, which are proportional to some power of a s (/z). 
Since the dominant part of the radiative corrections comes from a short distance region in the lattice four-momentum 
integral, we assume the scale /i to be 1/a. Secondly, discretization effects in the lagrangians will be of the order of 
(a|p|)™, where n is an integer and p is some soft momentum scale that characterizes the spatial momentum of the 
system. We take these soft modes to be of the order of Aqc d ■ Finally, there are power corrections to the heavy quark 
effective Hamiltonian, which are of the order of some power of Aqcd/M. 

In detail we expect the following corrections in our case, (i) Gluon and light quark actions: For the RG-improved 
gauge action, the leading discretization error is of the same order as the plaquette gauge action, which is (ciAqcd) 2 ■ 
For the 0(a)-improved Wilson quark action for light quarks, the leading error is of 0{a s aAqcD) since the coefficient 
csw is tuned at one-loop level only by the mean field improvement, (ii) Heavy quark action: For the 0(a)-improved 
Wilson quark action, the leading error appears in the 1/M 2 term at the tree level in the effective Hamiltonain, which 
is a source of systematic error of 0((Aqcd/M) 2 ). An additional error comes from the radiative correction that 
changes the relation between and Mb, and yields an uncertainty of 0(o. s Aqcd/M). (iii) Current corrections: 
The renormalization coefficient Za is computed only to one-loop accuracy, hence higher order uncertainties are of the 
order of a 2 . Other corrections to the current are present, but these are of the same order as those in the heavy quark 
effective Hamiltonian. 

The size of these corrections are estimated in Tables XIII and XIV for the B and D mesons. Numerical values are 
evaluated adopting the MS coupling at the scale fx = 1/a defined by (|30|) and Aqcd — 300 MeV, and substituting 
in M the physical B or D meson mass. In the case of Nf=Q, we only choose three representative (3 values as the 
variation of the errors is so mild across the available range of lattice spacings. The total uncertainty is estimated by 
adding all individual sources in quadrature. 

In Figure ^2] we replot fs d for Nf = 2 as a function of lattice spacing. The statistical error is shown by thick 
error bars, while thin lines represent the total error for which the statistical and estimated systematic errors are 
added in quadrature. The systematic error is the smallest at the finest lattice spacing ((3 — 2.1) as one can see in 
Table XIII. This confirms our expectation that the result from this (3 value provides the best estimate for fg d in the 
continuum limit. It is also important to note that the the data at coarser lattice spacings, especially that at (3=1.95, 
are consistent with the result at (3=2.1, if we take the systematic error into account. This bolsters our confidence in 
the estimate of the systematic error. 
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E. Continuum estimate 



In Figures |16(-|19| we plot our final results, including the estimated total error, for the continuum value of the 
heavy-light decay constants at a — 0. For the Nf = 2 calculation with dynamical quarks, the central value is taken 
from the data at the finest lattice spacing ((3—2.1), and the total error shown is obtained by quadratically adding the 
statistical and systematic errors. Numerically, we find for Nt — 2, 

f^r 2 = 208(10) (11) MeV, (33) 
f^= 2 = 250(10) (13) (^) MeV, (34) 

jBs ' = 1.203(29)(43)(+f), (35) 



JBd 

f2 =2 = 225(14) (14) MeV, (36) 
f%r 2 = 267(13) (17) (15°) MeV, (37) 

, \N f =2 

^ = 1.182(39)(46)(^ 1 ). (38) 

JDdJ 

The first error is statistical, and the second error is the cumulative systematic error outlined above. For the ratios, as 
ambiguities due to the renormalization coefficient arc eliminated, only the effect of the gluonic and light quark errors 
are included. In the case of those quantities involving the strange quark, the central value was taken from the strange 
quark mass defined from mx, while a systematic error was estimated from mass of the 4>. If instead of adding the 
systematic errors quadratically, we added them linearly, the final results, taking f Bd as an example, would be f Bd = 
208(10)(19) MeV. It is encouraging that our prediction for fn s with Nf — 2 is consistent with the recent experiments 
285±20±40 MeV (ALEPH §) and 280±19±44 MeV (CLEO ||). 

We also quote the results for the quenched case Nf = 0, for which we employ a constant fit to the five data points 
in the region a< 0.8 GeV -1 corresponding to (3 = 2.575-2.416. The estimated systematic error varies only slightly in 
this region, and we find for Nf = 0, 

f%r° = 188(3)(9) MeV, (39) 

c N f=° _ onr\(n\r-i K\(+6\ 



fB. 
fBd 



Id. 

fDd 



f^r = 220(2) (15) (I-) MeV, (40) 

N f =0 

= 1.148(8)(46)(±f), (41) 



/d<P° = 218(2) (15) MeV, (42) 
250(1) (18) (lg) MeV, (43) 

1.138(5)(45)(±g 9 ), (44) 



Da 

N f =0 



(45) 



where the systematic errors are assigned with the same strategy as for the case of Nf — 2. 

These quenched decay constants lie at the upper end when compared with those from previous quenched lattice 
calculations, whose recent summary is f Bd = 170(20) MeV, f B , = 195(20) MeV j§, and f Dd = 200(20) MeV, f D , 
= 220(^20) MeV J|]. One of the differences in this work from the previous ones is in the choice of the gauge action, 
with which the pcrturbativc matching factor is significantly different. Hence two-loop uncertainties of order 5% or 
so should be allowed for in the comparison. Another difference is the treatment of the 0(1/ M) correction arising 
from the field rotation (|l3|), which is included in our calculation but not in the previous studies. As we discussed in 



Sec. V A this term increases the value of decay constants. If we eliminate the 1/M correction in our calculation, we 
find /^J = ° = 182(3)(9) MeV = 213(2)(14) MeV, f^ f d =Q = 208(2)(14) MeV and /^ /=0 = 238(1)(17) MeV, i.e., 

a 3% effect for B and a 6% effect for D. Combined with uncertainties from the matching factors dicussed above, we 
consider that our quenched results with the RG-improved action are consistent with the previous data obtained with 
the plaquette gauge action. 
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F. Quenching effects 



/Bd 



N,=2 



In order to elucidate the effect of introducing sea quarks it is instructive to take the ratio of the results for Nt = 2 
and Nf = 0, for which we find : 

f N f =2 

(46) 
(47) 
(48) 
(49) 
(50) 
(51) 



fD. 



N f =2 




The errors quoted above are statistical only since we expect the systematic error to largely cancel in the ratio. We 
observe that the magnitude increases by 10-15% for the B meson decay constants when two flavours of dynamical 
quarks are introduced, which has statistical significance of 2 to 3 standard deviations. For the D meson decay constant, 
on the other hand, the observed increase is only 3-7%, and the effect is statistically less significant. For the ratio of 
decay constants we find only a small change from Nf — to Nf = 2. 

An increase of B meson decay constants in the presence of dynamical sea quarks has already been suggested in 
Refs. [ffiuj. Our results also show this trend, providing further evidence that the upward shift is real. 



VI. CONCLUSIONS 



In this paper we have presented a calculation of the heavy- light axial decay constants fsa, fBs, fDd, Ids and their 
ratios in lattice QCD with two degenerate flavours of sea quark (Nf = 2) where the same discretization scheme has 
been employed for the sea and light valence quarks. In order to carry out the calculation with the computational 
resources available, the heavy quarks are treated using an effective field theory approach, and the light quark and 
gluon fields actions are improved to minimize the discretization error. The calculation is also made in the quenched 
(Nf = 0) case since these decay constants have not been calculated before with this combination of actions. 

Our quenched continuum estimates for the decay constants are somewhat large compared with the current summary 
of quenched results. We note, however, that we have included an extra 1/M correction to the current, which was 
considered small and hence ignored in the previous calculations. We have demonstrated that this is not entirely true, 
giving a contribution of the order of +3% to fsd and /_b 3 and +6% to fud and fr>„. Allowing for the additional 
uncertainties from the use of one-loop perturbative renormalization factors, we consider our quenched results to be 
consistent with the previous results. 

In comparing our Nf — and Nf — 2 results we see that fsd and Jb s for Nf — 2 are significantly larger than the 
Nf = results by 2-3 standard deviations, indicating a shift of 10-15%. On the other hand, the same cannot be said 
for the decay constants fr>d and fo 3 - It is encouraging that our prediction fp a = 267(1 21) MeV with Nf — 2, where 
the total error is obtained by quadrature, is consistent with recent experiments. In conjunction with the available 
experimental data, our values for the Nf = 2 B meson decay constants f Bd = 208(15) MeV and f Bs = 250(±^|) MeV 
are consistent with the hypothesis that the Wolfenstein parameter p |}3| is positive. Given our results for Nf = and 
Nf = 2, it is reasonable to think that additional flavours of sea quarks will increase fg d and fs s still further, which 
in turn favours a positive value for p even more. 

The unsatisfatory aspect of our results is a sizable variation of the decay constants with lattice spacing. A possible 
origin of this problem is a necessity to include O(a) and higher improvement terms in the axial vector current. Higher 
order corrections in the renormalization constants may also be important at the coarse lattice spacings of a -1 1- 
2 GeV explored in the present simulation. The study of these issues is clearly needed to consolidate the results for 
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Nf — 2 and further explore the final goal of predicting the heavy-light decay constants for the realistic spectrum of 
dynamical sea quarks. 
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FIG. 1. Typical effective mass plots at /3 = f .8 for Nf = 2. The fit range is from 3 to 11. 



14 




15 




FIG. 3. Typical effective mass plots at /3 = 2.1 for Nf = 2. The fit range is from 5 to 21. 
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FIG. 6. The renormalization constant Za as a function of amo for the Wilson (at /3=5.9) and RG (at ,0=2.528) gauge actions. 
In the case of the RG action, Za is computed to specifically include the 1/M correction to the current while in the Wilson 
action it is not. The inverse lattice spacing is roughly 1.8 GeV (in the quenched approximation) for both cases. 
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FIG. 7. Ratio f^ otated / f^ nrotated _ i Q f decay constant including the current rotation to the lowest order current to that 
without the correction for Nf = 2 (circles and crosses) and Nf = (triangles). The gauge couplings were picked so that the 
lattice spacing roughly matched with each other. 
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FIG. 8. A plot of $(oM) vs 1/aM for Nf = 2 at /3=1.8. The data for four different sea quark masses are shown, and the 
light valence quark mass is set equal to the sea quark mass. 
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FIG. 9. Same as Fig. §, but for (3 = 1.95. 
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FIG. 10. Same as Fig. @, but for /3 = 2.1. 
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FIG. 11. A typical plot of $(aM) vs 1/aM for N f = at /3=2.334. 
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FIG. 14. A plot of $ s vs am q for (3 = 2.10 and AT/ = 2. 
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FIG. 16. fsd for Nf = 2 (filled circles) and Nf — (open circles) as a function of lattice spacing a. The error bar for the 
data points represents the statistical errors only, while those in the continuum limit (a—0) are the systematic and statistical 
errors added in quadrature. 
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FIG. 17. Same as Figure fifl but for f Be (K). 
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FIG. 18. Same as Figure [u], but for fn d ■ 
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FIG. 19. Same as Figure [T|, but for f Ds (K). 
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FIG. 20. A comparison of the ratio fBs(K)/fsd for Nf = and iV/ = 2. The error bars of the continuum limit results are 
the systematic and statistical errors added in quadrature. 
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FIG. 21. Same as Figure but for f D ,(K)/f D d- 
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FIG. 22. fs d with combined statistical and systematic errors for Nf = 2. The statistical error is shown by thick error bars, 
while thin lines represent the total error for which the statistical and estimated systematic errors are added in quadrature. 
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TABLE I. Simulation parameters for Nf 


= 2 lattices used 


in the heavy quark calculation. 


For the number of trajectories 


those in parentheses 


show the full ensemble 


generated. 









1.8 


3.162 


1.090(11) 


1.95 


2.812 


0.7882(85) 


2.1 


2.562 


0.559(11) 



TABLE II. Chirally extapolated parameters for Nf — 2 lattices used in the heavy quark calculation. The lattice spacing is 
fixed by p meson mass. 
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2.528 


1.355 


2.349 


0.612(6) 


195 


8 


2 


7.446(2) 


7.471(4) 


2.575 


1.345 


2.298 


0.574(6) 


200 


8 


3 


7.458(2) 


7.480(4) 



TABLE III. Simulation parameters for Nf = 0. The lattice size employed is 16 3 x 32 for /3=2.187-2.281 and 24 3 x 48 for 
f3=2. 416-2. 575. The lattice spacing is fixed by p meson mass. 
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p 


AT 


j- 

' * sea 


A 


D 
D 


1.8 


2 


0.1409 


1.09 


0.91 


1 o 
1.8 


o 
z 


n i /t on 
0.14JU 


i no 
1.09 


n ni 
0.91 


1 8 
1.8 


o 

z 


U. 1440 


i no 
i .uy 


u.y i 


l.o 


Z 


n i a Qn 
U.14oU 


i no 


n ni 
U.91 


1.95 


2 


0.1375 


1.28 


0.58 


1.95 


2 


n i on 
0.1o9 


1.28 


n k o 
0.08 


1.95 


2 


0.140 


1.23 


0.71 


1.95 


2 


0.141 


1.27 


0.6 


2.1 


2 


0.1357 


1.28 


0.54 


2.1 


2 


0.1367 


1.28 


0.54 


2.1 


2 


0.1374 


1.28 


0.54 


2.1 


2 


0.1382 


1.28 


0.54 


2.187 





_ 


1.28 


0.58 


2.214 







1.28 


0.58 


2.247 







1.28 


0.58 


2.281 







1.28 


0.58 


2.416 







1.28 


0.54 


2.456 







1.28 


0.54 


2.487 







1.28 


0.54 


2.528 







1.28 


0.54 


2.575 







1.28 


0.54 



TABLE IV. Smearing parameters used. The gauge-fixed smearing function takes the form Aexp (—Br) 
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p 


A 


Ai 


Bo 


Bi 


Co 


X 2 /d.o.f. 


2.187 


0.670(27) 


1.10(6) 


-0.79(8) 


-1.17(9) 


0.485(60) 


1.9/9 


2.214 


0.597(29) 


1.01(9) 


-0.66(8) 


-0.95(13) 


0.373(55) 


1.3/9 


2.247 


0.556(25) 


1.00(8) 


-0.57(7) 


-0.88(11) 


0.296(51) 


0.7/9 


2.281 


0.480(25) 


1.00(13) 


-0.43(6) 


-0.86(17) 


0.196(36) 


0.9/9 


2.334 


0.412(22) 


0.93(10) 


-0.36(5) 


-0.78(14) 


0.173(31) 


2.1/16 


2.416 


0.319(12) 


0.84(5) 


-0.20(2) 


-0.60(5) 


0.055(10) 


0.4/11 


2.456 


0.317(2) 


0.52(24) 


-0.21(3) 


-0.31(22) 


0.065(7) 


0.7/11 


2.487 


0.264(10) 


0.78(4) 


-0.16(2) 


-0.47(4) 


0.044(6) 


0.3/11 


2.528 


0.237(16) 


0.82(34) 


-0.13(1) 


-0.50(23) 


0.032(4) 


1.2/11 


2.575 


0.235(9) 


0.67(5) 


-0.14(1) 


-0.36(4) 


0.043(4) 


3.8/11 



TABLE V. Chiral HQET fit parameters for N f = 0. 



p 


A 


Bo 


Co 


X 2 /d.o.f. 


2.187 


0.785(23) 


-0.89(7) 


0.456(63) 


0.7/4 


2.214 


0.699(22) 


-0.75(7) 


0.366(57) 


0.6/4 


2.247 


0.650(20) 


-0.65(6) 


0.293(51) 


0.3/4 


2.281 


0.583(15) 


-0.54(5) 


0.228(38) 


0.8/4 


2.334 


0.492(13) 


-0.42(4) 


0.172(32) 


0.4/4 


2.416 


0.365(15) 


-0.24(3) 


0.063(22) 


0.006/4 


2.456 


0.357(10) 


-0.25(2) 


0.076(11) 


0.1/4 


2.487 


0.315(8) 


-0.19(2) 


0.045(8) 


0.03/4 


2.528 


0.288(12) 


-0.16(2) 


0.030(13) 


0.2/3 


2.575 


0.273(7) 


-0.17(1) 


0.051(6) 


0.4/3 



TABLE VI. Strange (defined from the K) HQET fit parameters for N f = 0. 



p 


A 


Bo 


Co 


X 2 /d.o.f. 


2.187 


0.814(22) 


-0.95(7) 


0.492(60) 


0.7/4 


2.214 


0.723(20) 


-0.80(6) 


0.392(55) 


0.7/4 


2.247 


0.667(20) 


-0.67(6) 


0.300(52) 


0.3/4 


2.281 


0.611(15) 


-0.59(5) 


0.257(46) 


1.1/4 


2.334 


0.512(13) 


-0.45(4) 


0.188(28) 


0.5/4 


2.416 


0.380(13) 


-0.25(3) 


0.063(18) 


0.006/4 


2.456 


0.362(7) 


-0.25(2) 


0.076(10) 


0.1/4 


2.487 


0.326(8) 


-0.20(1) 


0.047(8) 


0.02/4 


2.528 


0.301(20) 


-0.17(4) 


0.036(22) 


0.07/4 


2.575 


0.282(6) 


-0.18(1) 


0.052(6) 


0.3/4 




TABLE VII. 


Strange (defined from the 4> 


HQET fit parameters for N f = 0. 
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p 


A 


A 1 


A 2 


So 


Si 


Co 


X 2 /d.o.f. 


1.80 


0.95(2) 


0.48(4) 


-0.10(2) 


-1.27(6) 


-0.31(3) 


0.76(3) 


30.4/26 


1.95 


0.51(2) 


0.61(6) 


-0.14(4) 


-0.72(8) 


-0.49(7) 


0.58(8) 


18.7/26 


2.10 


0.25(1) 


0.56(9) 


-0.34(15) 


-0.18(1) 


-0.25(2) 


0.06(5) 


7.8/26 



TABLE VIII. Chiral HQET fit parameters for N f = 2. 



a 


A 


Ai 


A 2 


Bo 


Bi 


Co 


X 2 /d.o.f. 


1.80 


1.09(2) 


2.57(28) 


-6.83(119) 


-1.43(4) 


-1.14(18) 


0.78(3) 


33.4/26 


1.95 


0.60(2) 


2.58(37) 


-4.71(225) 


-0.79(7) 


-2.06(43) 


0.57(6) 


19.8/26 


2.10 


0.30(1) 


1.32(36) 


-3.86(287) 


-0.20(1) 


-0.58(10) 


0.06(5) 


6.2/26 



TABLE IX. Strange (defined from the K) HQET fit parameters for Nf = 2. 



p 


A 


Ai 


A 2 


Bo 


Bi 


Co 


X 2 /d.o.f. 


1.80 


1.13(2) 


2.64(26) 


-7.25(156) 


-1.49(3) 


-1.13(18) 


0.81(3) 


30.9/26 


1.95 


0.62(3) 


2.65(37) 


-5.55(217) 


-0.82(8) 


-1.97(42) 


0.59(7) 


11.2/26 


2.10 


0.31(1) 


1.22(36) 


-2.69(282) 


-0.21(1) 


-0.63(9) 


0.07(4) 


4.5/26 



TABLE X. Strange (defined from the <j>) HQET fit parameters for N f = 2. 



n f 


p 


f Bd (GeV) 


f B .{K) (GeV) 


f B .{4>) (GeV) 


fBs(K)/f Bd 


fB,{4>)l 'fBd 


2 


1.80 


0.287(7) 


0.331(5) 


0.340(5) 


1.152(19) 


1.181(26) 


2 


1.95 


0.234(8) 


0.276(7) 


0.283(8) 


1.179(43) 


1.211(45) 


2 


2.10 


0.208(10) 


0.250(10) 


0.258(10) 


1.203(29) 


1.241(36) 





2.187 


0.229(7) 


0.268(5) 


0.276(5) 


1.171(13) 


1.121(16) 





2.214 


0.220(8) 


0.258(5) 


0.265(5) 


1.169(20) 


1.202(25) 





2.247 


0.223(6) 


0.260(5) 


0.266(5) 


1.165(16) 


1.194(19) 





2.281 


0.204(7) 


0.244(4) 


0.254(3) 


1.196(31) 


1.243(40) 





2.334 


0.195(7) 


0.232(4) 


0.240(4) 


1.186(29) 


1.227(35) 





2.416 


0.188(6) 


0.212(5) 


0.220(5) 


1.126(11) 


1.169(15) 





2.456 


0.204(11) 


0.228(4) 


0.232(3) 


1.111(58) 


1.103(39) 





2.487 


0.184(5) 


0.217(4) 


0.226(3) 


1.182(14) 


1.221(18) 





2.528 


0.182(11) 


0.219(6) 


0.228(9) 


1.190(104) 


1.235(130) 





2.575 


0.192(6) 


0.221(4) 


0.227(4) 


1.151(18) 


1.183(21) 






TABLE XI. 


Decay constants fsd and /b, 


at each bare gauj 


;e coupling. 
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n f 


p 


f Dd (GeV) 


fos(K) (GeV) 


f D .{4>) (GeV) 


f D s{K)/f Dd 


fnsW/fn, 


2 


1.80 


0.301(9) 


0.346(7) 


0.352(6) 


1.151(18) 


1.168(24) 


2 


1.95 


0.284(13) 


0.312(11) 


0.316(11) 


1.097(40) 


1.111(39) 


2 


2.10 


0.225(14) 


0.267(13) 


0.277(13) 


1.182(39) 


1.223(47) 





2.187 


0.258(4) 


0.295(3) 


0.300(3) 


1.143(7) 


1.165(8) 





2.214 


0.247(4) 


0.284(3) 


0.290(3) 


1.149(9) 


1.172(12) 





2.247 


0.250(4) 


0.287(3) 


0.293(3) 


1.146(8) 


1.169(10) 





2.281 


0.236(4) 


0.274(2) 


0.282(2) 


1.160(17) 


1.191(22) 





2.334 


0.229(4) 


0.263(3) 


0.269(3) 


1.148(12) 


1.175(14) 


n 
u 


2.416 


n 9991' 41 








1 1 1 1 1 





2.456 


0.229(8) 


0.254(3) 


0.257(3) 


1.129(73) 


1.119(46) 





2.487 


0.215(3) 


0.249(2) 


0.256(2) 


1.159(7) 


1.188(9) 





2.528 


0.215(10) 


0.250(4) 


0.256(6) 


1.133(67) 


1.159(81) 





2.575 


0.215(4) 


0.245(3) 


0.251(3) 


1.138(13) 


1.163(16) 






TABLE XII. 


Decay constants /dj 


and /ds at each bare g 


;auge coupling. 





P 


1.8 


1.95 


2.1 


{Kqcd/MbY, (Aqcd/MdY 

2 

(ciAqcd) 2 

CtsClAQCD 

oi s Kqcd/Mb, ol s Kqcd /M d 


<1%, 3% 
5% 
11% 
6% 
1%, 3% 


<1%, 3% 
3% 
6% 
4% 
1%, 3% 


<1%, 3% 
3% 
3% 
2% 
1%,3% 


total (linear) 
total (quadratic) 


23%, 28% 
14%, 14% 


14%, 19% 
8%, 9% 


9%, 14% 
5%, 6% 


TABLE XIII. An estimate of systematic 


errors for Nf —2. 


Aqcd is taken to be 300 MeV. By convention, 


we use the 



quadratically summed error. 



p 


2.187 


2.416 


2.575 


(Aqcd/MbY, {Kqcd/MdY 

(ciAqcd) 2 
q s oAqcd 
oi s Aqcd /M b , o-sAqcd /M d 


<1%, 3% 
3% 
9% 
5% 
1%, 3% 


<1%, 3% 
3% 
5% 
3% 
1%, 3% 


<1%, 3% 
4% 
3% 
3% 
1%,3% 


total (linear) 
total (quadratic) 


18%, 23% 
11%, 12% 


12%, 17% 
7%, 8% 


11%, 16% 
5%, 7% 


TABLE XIV. An estimate of systematic 


errors for Nf =0. 


Aqcd is taken to be 300 MeV. By convention, 


we use the 



quadratically summed error. 
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